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0.  Introduction 

the  problem  of  selecting  populations  better  than  a  control  with 

respect  to  a  location  parameter  under  ordering  prior  has  been  considered 

in  [4].  In  this  paper  we  consider  the  case  of  binomial  populations 

(important  for  discrete  data)  for  which  the  parameters  of  interest  are 

not  the  location  parameters  as  was  the  case  studied  in  [4],  We  consider 

both  cases  when  the  parameter  of  control  is  known  and  unknown.  For  the 

case  of  known  control,  we  propose  an  isotonic  procedure  which  is  given 

in  Section  2.1.  The  results  in  this  section  deal  with  both  cases, 

namely,  the  sample  sizes  all  equal,  and  unequal.  Some  recursive 

relations  are  derived  for  computing  the  constants  required  for  the 

proposed  procedure.  When  the  control  is  unknown,  a  conditional 

isotonic  procedure  is  proposed  in  Section  2.2.  This  procedure  provides 

a  conservative  solution  for  the  unconditional  procedure.  Brief  tables  of 

associated  constants  for  the  proposed  procedures  are  given  in  Table  1  and 
Table  2A 

1.  Notations  and  Definitions 

Assume  Hq,^,...,^  are  all  binomial  populations  such  that  has 
density  b(m;  p^),  1  =  0,l,2,...,k.  It  is^assumed^Kat  p^  <  Pg  <...<_  p^, 
however,  the  actual  values  of  these  p^'s  are  not  known.  We 

♦this  research  was  supported  by  the  Office  of  Naval  Research  contract 
N00014-75-C-0455  at  Purdue  University.  Reproduction  in  whole  or  In 
part  is  permitted  for  any  purpose  of  the  United  States  Government. 


2 


consider  Hq  as  a  control  and  our  goal  is  to  select  a  subset  of  these 
k  populations  so  that  all  "good"  populations  are  included  in  the  subset 
selected,  where  is  considered  "good"  if  and  only  if  p^  >  Pq.  By  a 
correct  selection  (CS)  we  mean  the  selection  of  any  non-trivial  subset 
which  contains  all  good  populations. 

Let  { (  Pq.P-j  ,  • .  •  *P|^)  1 0  <  Pj  i  P2  — '  **.5.P|(  <  1 ®  <  Pq  <  11* 

Let  us  denote  the  sets  ai  =  {i,  i+1 . k},  1  <  i  <  k  and  aQ  =  4  (the 

empty  set).  If  action  ai  is  taken,  it  means  the  subset  »irj+j . 1^} 

of  the  k  populations  is  selected.  Since  by  our  assumption  that  p.  are 
ordered  according  to  an  ascending  (simple)  ordering  prior,  it  is 
therefore  appropriate  to  restrict  ourselves  to  the  action  space 

G  *  {aQ,a1 ,a2,. . . ,ak>.  For  given  positive  integer  ng.n^ . n^,  we 

assume  n^  independent  observations  are  drawn  from  ir^  (i  =  0,l,2,...,k) 
which  are,  respectively,  denoted  by  X11,X12,...,Xi  .  The  sample 

space  is  denoted  by 


Z  =  (X  e  I 


W* 


( ,x. 


01  ^02’ 


where  I  denote  the  set  of  non-negative  integers. 


We  restrict  ourselves  to  isotonic  selection  procedures  S 

which  satisfy  the  P*-condition,  i.e.  inf  P  (CS|fi)  >  P*  where  P*  is  a 

0€fi  - 

preassigned  value. 

A  poset  (S,  <)  denotes  a  non-empty  set  S  with  a  binary  partial 
order  <  defined  on  it.  A  real-valued  function  f  defined  on  a  poset 
(S,  <)  is  called  isotonic  if  f  preserves  the  order  on  S,  i.e.  x  <  y. 
Implies  f ( x)  £  f(y).  Let  g  be  a  real -valued  function  and  let  W  be  a 
positive- valued  function,  both  defined  on  a  poset  (S,  <).  An  isotonic 
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function  g*  on  S  Is  called  an  Isotonic  regression  of  g  with  weight  W  If 

l  [g(x)-g*(x)]2W(x)  attains  Its  minimum  values  over  set  of  all 
x€S 

Isotonic  functions  on  S.  It  Is  well-known  (see  [2])  that  there  exists 
one  and  only  one  Isotonic  regression  of  a  given  g  with  a  given  weight 
W  defined  on  S.  Some  algorithms  such  as  "pool-adjacent-violators"  or 
the  so-called  "up-and-down  blocks"  are  referred  to  In  [1]  and  [2]. 

Let  nQ  =  0  and  n-j  *  ng  =...=  nfc  =  n  and  m  *  1.  Let  5^  denote 
the  sample  mean  from  ir. ,  1  =  l,2,...,k.  The  Isotonic  regression  of 
X.  with  common  weight  n  is  a  maximum  likelihood  estimate  for 
P  =  (p1 ,p2«* • • .Pk)  which  is  given  by  the  following. 

Theorem  1.1  ([2]  p.  102} :  The  maximum  likelihood  estimate  for 

p  =  (pj ,P2«* • • ,Pk)  with  p1  <  Pg  <...<  pk  is  given  by  the  isotonic 

regression  of  5^  with  common  weight  n,  l.e.  X  =  (Xj  k,X2  k . Xk  k) 

k  _  o 

minimizes  7  ( x.-p . )  n  where,  by  the  max-mi n  formula  of  [1],  we 
i=l  1  1 

have 

(1.1)  X,  Li  =  max  {  min  ((X+X..,  +...+  X+)/(t-s+l))} 

1,K  l<s<i  s<t<k  s  s  1  1 


=  max  X.  l 
1 <s<1  S,K 


where 


0.2)  X$>k  ■  min{X$,  \  (Xs+Xs+1 (Xs+Xs+i  +...+  Xk)/(k-s+l ) > . 


2.1  Isotonic  Selection  Procedure  <5^  When  pQ  is  Known, 
a.  Equal  Sample  Size  Case 

Since  Pg  is  known,  we  take  ng  =  0;  assume  n^  *  n2  =...=  nk  3  n. 
Without  loss  of  generality,  we  may  consider  m«  1,  l.e.  Is  a 
Bernoulli  (p^).  For  given  positive  constants  dj,  d2,...,dk 
(0  <  d^  <  Pg,  to  be  determined  later),  we  propose  the  following 
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unconditional  procedure  <5^: 


«•, (X)  =  ag(X),  where 


€(X)  =  min{i|X.j  k  iPg-d^}  and  X^  k  is  defined  by  (1.1). 
Since  6-j ( X)  depends  on  vector  d  =  (d-j  ,dg». .. .d^),  we  may  denote 
it  by  <S -j (d)  when  there  is  no  confusion. 


Determination  of  d  for  6 ^ ( d) 

In  order  to  satisfy  the  basic  P*-condition,  we  need  to  compute 

inf  P(CSj61(d)).  For  notational  conveniences,  we  define 
n 

fl1  =  IP  €  Q|pk_i  <  pQ  <  Pk_i+i>.  1  =  1.2 . k-1 


nk  =  (p  €  n|pQ  <  p^ 


and 


n0  -  (P  e  fl|pk  <  p0>. 

k 

Then,  fi.  are  disjoint  and  n  =  U  fl. .  Again  for  notational  convenience, 
1  i=0  1 

^  a 

when  there  is  no  confusion,  we  denote,  respectively,  X.  t  and  X.  t  by 

•  I1'  J 

A  A 

X^  and  Xj  for  a  given  fixed  k.  Then,  for  any  p  E  s., 


k-i+1 


(2.1)  Pp(CS|61(d))  =  P?{  ^  tXj  iP0-dj}> 


=  P 


k-i+1  j  x 


{  U  U  (Xr  >  p0-dj» 
c  j*l  r=l  J 


-  Pp{xk-i+l  -  p0'dk-i+l} 

Since  pp(Xk.i+1  LP0-dk_i+i)  is  increasing  in  Pk_1+j  for  j  =  l,2,...,i, 
keeping  all  other  (i-1)  parameters  fixed,  hence,  the  right  hand  side 
of  (2.1)  attains  Its  minimum  at  Pk_^  =  Pk_i+2  =***=  Pk  =  Pq*  0n 
other  hand,  if  we  take  a  special  vector 


■t . , 
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k-i 

(2.2)  pQ  =  (p0,0,0 . O.Pq.Pq,. . . , Pq )  €  the  closure  of  ttj. 


then,  we  see  that 


P 


k-i+1 

Po(  j=l 


^(xr  >  P0'dj))  =  Pp0(Vl+l  -  p0“dk-i+l} 


since  X.  =  0  a.s.  when  p.  =  0  for  j  =  1,2,..., k-i.  Hence, 

J  j 

inf  P(CS|«,(d))  attains  at  p  =  pQ  as  defined  by  (2.2).  Again,  since 


(2.3)  inf  P  ( CS | 6 , ( d) )  =  min  {inf  Pn(CS|6,(d))> 

n  P  1  "  l<i<k  n.  P  1  " 

because  for  p  €  Jig,  any  action  in  a  is,  according  to  our  definitions,  a 
correct  selection.  Therefore,  if  for  each  i  (1  <_  i  <_  k)  inf  P  (CS|6,(d)  >_  P*, 

a  *  P 

then  the  P*-condition  holds  for  (d) .  Now, 

(2.4)  1nfP{CSH,).P  (XMtl  >  Vd,.w) 

ni  ~u 

where 

3q  -  tPg *Pq * •  •  *  »Pg)  ^  (0»1)  ,  and  oj  = 

It  follows  from  (1.2)  and  (2.4)  that 


(2.5) 


inf  P{CS|61 (d)>  =  Pr{Y]  > 
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■^(Yl+Yg)  >_  c^ , . . . , 

|(Y1+Y2+...+Yi)>  c^) 


where  Y^ ,  Y2,...,Y^  are  i.i.d.  with  Y^  being  b(n;  pQ)  and 
ci  "  n^Po“ai^*  Hence»  it  suffices  to  compute 


i  a 


...nil 
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(2.6)  a.(a.)  =  Pr { Y,  >  er  Y]  +  Y?  >  e2....,  j  Yj  >  e^ 

j  1 

where 

(2.7)  =  j  c.  =  jn(pQ-a.),  j  =  1,2,. ..,1. 

Define 

(2.8)  «j(«,)  -  Pr(Y,  i*,.  V,  *V2>e2 . J-1.2 . 1. 

Letting  <a>  =  -[-a],  I.e.  the  smallest  integer  no  less  than  a,  we  have 
the  following  useful  lemma. 


Lenina  2.1.  (1)  a^)  =  l  (J)pS(l"P0)n’r  where 

r=rQ 

r0  a  <n(P0-a1)> 

(ii)  a.(a.)  =  l  g(r.){  l  g(r2){  £  g(r3)  ... 

J  r1=<c>  r2=<2c-r1>  r3=<3c-r^-r2> 

{  l  g(rj},  j  =  1,2 . 1 

rja<^C"rl"r2* '  *rj-1> 

where 

c  *  n(pQ-a.j)  and 
9(r)  .  (;)p0r(l-p0)n-r. 


Proof:  To  compute  (1)  Is  straightforward.  To  prove  (1i),  define 

A^o.e)  =  Pr{Y,  >  a,  li  +  Y.  >  a  +  8....,  f  Y  >  a  +  (j-1)8l 

J  lie.  r=l 

where  YpY2,...,Yj  are  lid  b(n;  Pg).  Conditioning  on  Y^  =  r-j,  we  obtain 


Aj(«.f»)  =  l  (")Po<1“P0)  Pr*(Y2  ^  a+6_rl‘  Y2+Y3  i  a+2^r] 


r=<a> 


•  ••»  £  Vr  >  a  +  ( J-l )&-r-| > 


l  (>0<'-P0>""V  l(»+e-ri.S). 
=<«>  U  J 


r^=<a> 


Taking  a  =  3  =  n(pg-d)  and  using  (i)  and  by  mathematical  induction, 
we  obtain  (ii). 

Hence,  for  given  a.,  this  lemma  gives  a  direct  method  of  computing 
a^(a^).  For  some  special  values  of  n(=  5(1)10),  p0(=0.1(0.1)0.5) , 

P*(=  0.90,0.95)  and  i(=  1(1)4),  the  smallest  values  of  satisfying 
a. (0^)  >_  P*  are  tabulated  in  Table.  It  is  to  be  noted  that  dfc_1.+.  = 

(i  =  l,2,...,k)  in  procedure  6-j (d) . 

Now  define, 

aj(a,g)  =  P{n(pQ-a)  <  Y1  <  n(pQ-$),  2n(pQ-a)  <  Y-j+Yg  <  2n(pQ-B 
...,jn(p0-a)  <  |  Yr  <  jn(p0-3)}. 

Then,  from  analogous  arguments  as  in  Lemma  2.1,  we  have 

[c£]  C^Cg- ^ -j  D  C3c2~r  1  “rgD 

Corollary  2.1:  aUa.B)  =  l  g(r,){  f  g(r2){  £  gi 

J  rla<cl>  r2=<2c1-r1>  r3=<3c1-r1-r2> 


t  1  g(^)>> 


where 


rrroo 


C1  =  n(p0-o) 
c2  -  n(p0-3) 

r0j  =  <jcrr1-r2-...-rJ.1>,  j  =  2,3,... 
r0J  *  tjc,-rrr.-...-r,  ,]. 
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b.  When  the  sizes  are  not  necessarily  equal.  We  also  take  ng  =  0 
since  Pq  is  known,  and  assume  m  =  1.  Then  the  isotonic  estimates 
in  (1.1)  and  (1.2)  become 

X.  k  =  max  X  . 

’  l<s<i  S,K 

where 

(2-9)  *s,k  =  min{*s*  ("SXS  +  ns+l\+l^ns+ns+l^  *' '  ” 

(nsWl^l+-  *  *+\V^ns+ns+l+*,*+nk^' 

For  our  notational  simplicity,  we  define 
i  =  k-i+1  and 

(2.10) 

mi,j  =  ni+ni+l  +“-+  ni+j-l‘ 

Then  (2.5)  becomes 


(2.11)  inf  P  { CS 1 6 1 }  =  P  {Z1  Cj,  Z^+Z2  c2,. . .  ,Z^+Z2+. .  .+Z.  >  c.) 

where 

C1  =  mi,i(Po-ai) 

c2  =  mi  ,2^p0-ai  ^ 

(2.12)  . 

ci  =  mi,i^0-ai) 

and  Zr  Z2,...,Zi  are  iid  with  Z.  being  b(nk_i+j;  pQ).  Then,  (2.11) 
can  be  computed  according  to  the  following.  Define 

(2.13)  bj  |(c)  =  ,  Z}+12  c2»...»  ^ j ^ ^  ~  1 ,2, . . .  ,i , 

where  c  *  (c1  ,c2,. . .  ,ci )  and  Z^s  are  defined  by  (2.12).  Then,  we  have 


=  P  L  >9(Vl*Trl){r  L  r  g(nk-i+2‘r2>‘ 

rj! 

....  i  9(nk_i+j »rj) J  j  -  1,2,... ,i; 
rj=rjO 

where 

9(n;  r)  =  -P0)n_r 

rj0  =  <cj-rrr2"-'*‘rj-l> 

rjl  =  nk-i+j 

and  c^Cg,...,^.  are  defined  by  (2.12)  and  (2.10). 

Proof:  Conditioning  on  21  =  r-j  in  (2. 13)  and  following  analogous 
arguments  as  those  in  Lemma  2.1,  we  obtain  the  result. 

Corollary  2,2.  If  for  given  P*,  the  constants  d  =  (d^.dg,.. .,dk) 

associated  with  6,  ( d)  are  so  chosen  that  b-  .-(c)  >  P*  for  all 
'  ~  *  *  *  - 

i  =  l,2,...,k,  then  pp( CS 1 6-j (d)  >  P*  Vp  6  n,  where  c  and  b.  ^ ( c) 

are  defined  by  (2.12)  and  (2.13),  and  dk_i+1  =  ,  i  =  1,2 . k. 

Proof:  The  result  follows  from  (2.3),  (2.11)  and  (2.13). 

Computations  of  a.  in  Table  1 

For  given  n,  i,  pQ  and  P*,  we  start  by  taking  v1  =  pQ  -  Usin 
Lemma  2.1,  we  can  compute  a^v^.if  a^v^)  <  P*,  we  take 

"i  =  vl  +  TH  =  Po»  otherwise,  we  take  v2  «  V]  -  ^  =  pQ  -  And 

compute  ai(Vg).  If  a^Vg)  <  P*.  we  take  a.  =  v?  +  i  =  pQ  - 

otherwise  take  v^  =  v2  -  This  process  continues  until  for  the 

first  time  a^(vr_^)  >_  P*  and  ( vf)  <  P*.  Then,  we  stop  and  take 
r—  1 

=  v  .  =  p„  -  -rr.  The  reason  that  each  time  we  decrease  the  amount 
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of  jjj-  is  that  (2.5)  remains  unchanged  as  long  as  the  value  of  ici 
changes  by  an  amount  less  than  1,  i.e.  in^Q-o^. )  <  1. 

For  computations  of  a  and  B  such  that  a^a.B)  >  P*  given  by 
Corollary  2.1  are  analogous. 


2.2.  Selection  Procedure  ^  and  Conditional  Selection  Procedure  6^ 
When  Pq  is  Unknown. 


For  simplicity,  we  assume  m  =  1  and  nQ  =  n1  =...=  n^  =  n.  When 
P0  is  unknown,  we  propose  the  following  procedures. 

52<?;  !>  •  ^(x;  u)  "here 

e^x;  u)  =  min{i(Xij(c  >  Xq-u.^ 
where 

X-  k  is  defined  by  (1.1)  and  (1.2). 

We  also  propose  a  conditional  procedure  63  as  follows. 


s3<?;  ?l?>  *  ac2(x;  v|t) 
€2(x;  v 1 1)  =  min{i  |X..  ^  > 


where 


Vvi 


i 

given  that  J  X.  =  t • } 

j=0  J  1 


where 

j 

t  -  ( t -j ,  1 2 » •  • . » t^ ) ,  tj  —  ^  X^.,  ^  -  ^»2».**>k. 

We  note  that  when  (Xq,X-j ,. . . ,X^)  is  observed,  t  is  automatical ly 
known.  However,  in  some  situations,  the  experimenter  may  have  values 
of  X.  ^  -  Xq  without  knowing  (XQ,X^ ,.. . ,X^).  For  example,  the  given  data 
might  be  X-|  -  Xq,  X2  -  Xq,...,X^  —  Xq. 

For  our  convenience,  we  denote  62  ant*  respectively,  by  62(u) 
and  «_ ( v I t )  when  there  is  no  confusion. 

! 

t 

l 


i 
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For  fixed  i,  and  any  p  € 

k-i+1  j 


Pn(CS |62(u) )  =  P  (  U  U  (X  >  -  u.)) 

P  t  -  ?  j=l  r=l  r  u  J 

-  Pp{*k-i+l  -  Vuk-i+l}* 


For  fixed  pQ,  the  right  hand  side  is  an  increasing  function  of  Pk_1+1 
keeping  all  p.  .  .  (2  ^  j  <_  i)  fixed.  By  the  same  arguments  as  in  the 

K  “  I  *  J 

last  section,  we  see  that  the  right  hand  side  attains  its  minimum  at 
pk-i  +  l  =  pk-i+2  ="‘=  pk  =  p0*  By  choosing  Po  defined  bY  (2-2)»  and 
applying  the  analogous  arguments,  we  conclude  that 

(2.14)  inf  P  (CS|«2(u))=  inf  P{Y1  >  Yg-Wj,  Y1+Y2  >  2(Y0-w2), 


°i 


OiPqI1 


...,Y1+Y2+...+Yi  >1(Y0-W.)}. 

where  Yg,Y,  ,Y2>. . . ,Y..  are  iid  such  that  YQ  is  b(n;  pQ)  and 


Wj  =  j  n  u.. 


(2.15) 

Now  define 

(2.16)  A(i ;  ti,ui)  =  (Y]  >  Yq-w,,  Y1+Y2>2(Yq-w2) . Y^Y^.+Y.  ^(Yq^J, 

1  Yj  -  V 

j=0  J 

(2.17)  B(i ;  u.)  =  {Y,  >  Yg-Wj ,  Y1+Y2>2(YQ-w2), . . .  .Y^Y^. .  .+Y.  >  KYg-w^) 

where  wi  is  defined  by  (2.15).  Then,  we  have 

n(i+l)  i  i 

(2.18)  P(B(i;  u. ))  =  l  P(B(i;  u) |  l  Y.  =  t  )P(£  Yj  -  t,) 

1  t.=0  j=0  J  j=0  J 

n(i+l)  P(A(i ;  t.,u))  i 

=  t?0  "(i+1>  nl+n-t.  P(  iQ  Yj=ti}- 

'i0  (  t.  >P0  (H>0)  1 

If  for  any  0  iPgl  and  any  t^  (0  _<  t^  _<  (i+1  )n) 

n(i+l)  t.  ni+n-t. 

(2.19)  P(A(i;  ti,ui))>  P*(  t  )Pg  O-Pg) 


holds,  then  it  follows  from  (2.14),  (2.18)  and  (2.19)  that 
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(2.20) 


Define 


inf  P  (CS|«2(u))  >  P*. 


(2.21)  ^(t.s.u)  =  t(x1,x2,....xi)|x1  i  s-nu,  x^+x2  >  2s-2nu, 

i-1  1  i 

J.  Xj  1  (i-Ds-(i-l)nu^  X.=t>  c  I 


(2.22)  q^t.s.u)  «  I  ("  >(£ 

1  ^-(t.s.u)  1  2  i 

(2.23)  (t ,u)  =  {(x0,x1,x2,...,x.)|x1  1  ynu,  X]+x2  >  2xQ-2nu 

-j,xj  i  “o"'™' 

(2.24)  9l(t.U)  -  1  ("  )("x  >•••(")• 

1  ^(t.u)  0  1  1 


Then,  we  have 

t.  ni+n-t.  , ni+n, 

(2.25)  P(A(i;  t.  ,u))/(nJ_n)p()1  (1-Pq)  =  9,(^1  ,u)/(  t_  ). 

It  follows  from  (2.19)  and  (2.25)  that  in  order  to  find  ut  so  that 
62(u)  satisfies  the  P*-condition,  it  suffices  to  find  such  that 

for  all  ti  (0  <  ti  <  n(i+l ) ) 

ni+n 

(2.26)  g^t^,  u^  >  P*(  t  ) 


holds. 


To  compute  gi (t ,u)  for  given  t  and  u,  we  may  apply  the 


following  theorem: 


r 


Define 


1-1 

(2.27)  ^(t.a.fl)  =  {(xrx2 . xi)|x1  >  a,  Xj+Xg  >  a+0 . ^  *j 

i  i 

>  a+(1-2)0,  ^  Xj  =  t>  c  I  . 

(2.28)  (t  ,y)  »  {(x0,x1,...,x.)|x1  >  x0-y,  x^x2  >  2(xq-y)  , . . . , 

j  xj  ±'(VxK  l  xj  =  11 


(2.29)  u,(t.«.9)  -  l  00-0 

€i(t,a,8)  1  2  *1 


(2.30)  v . ( t ,y)  -  i  <;  x;  >...<;  >• 

1  ?i(t,Y)  X0  xl  Xi 

Then,  we  have  the  following 

n  n 

Theorem  2.2:  (i)  u,-(t,a,0)  =  £  (r  H  [ 

-  1  r-j=<a>  1  r2=<a+0-r-,>  2 

(  I  (J.)...*  I  <"  H.-r.-r  " 

rj’xj  3  ri-l  xi-l  '  2  11 


where 


Xj  =  <a+(j-l)0-rj-r2-. . . - r j _ -j >  j  =  2,3,...,i-l, 

X  =  min{n,t-r^-r2-.  ..-r^g} 

S2 

(ii)  v. ( t iy)  =  l  (")u,(t-s,a,s-Y)  "here 

s-s, 

=  max{0,t-in) 
s2  =  m1n(n,t, 

and  u^t-s.a.s-Y)  is  given  by  (i). 


Proof:  To  show  (i),  we  note  that  by  conditioning  on  x1  =  r^  in  the 

set  £.(t,a,8),  and  the  lower  bound  of  r-j  is  at  least  a,  we  thus 

n 

obtain  u^t.o.e)  ■  l  ("  )u.  ,(t-r.,  a+B-r,,s).  Using 

rr<a>  1  "  r0 

mathematical  induction  and  noting  that  u2(t,a,8)  =  £  (J!)(J]r)  where 

r=<a> 

rQ  should  exceed  neither  n  nor  t-s,  we  thus  obtain  (i). 

To  show  (ii),  we  condition  on  the  value  of  Xq  by  taking  Xq  =  s  in 

the  set  c_.(t,y).  Then,  in  the  set  ^(t.y),  we  have  simultaneously  the 
i  i 

conditions  |  >  i(s-y)  and  ^  xj  =  t-s*  In  orcler  that  ^(t.y)  is 

non-empty,  we  should  have  t-s  _>  i(s-y)  and  this  determines  the  upper 

bound  of  s.  On  the  other  hand,  we  also  note  that  s  should  exceed  neither 

n  or  t  and  this  concludes  the  value  s2  in  (ii).  For  the  lower  bound,  we 

note  that  s  should  not  be  less  than  t-in  if  this  is  positive.  Finally, 

i 

we  note  that  the  two  conditions  Ex.  >  i(s-y)  and  I  x.  =  t-s  combine  into 

J  |  J 

i 

|  Xj  =  t-s  and  it  concludes  (ii). 

This  completes  the  proof  of  the  theorem. 


Now,  by  taking  a  =  s-nu  and  b  =  s-nu,  (t,s-nu,s-nu)  =  <(>.(t,s,u), 
also,  taking  y  *  nu  we  have  5^(t,nu)  =  ^(t.u).  This  concludes  the 
followi ng 

Corollary  2.3:  (i)  q.(t,s,u)  =  l  (J  ){  l  ("  ){... 

r^=<a>  1  r2=<2a-r.j>  2 


where 


<  i  <".)--•{  I  <"  >i 


rj'xj rj 


t-r,-r2-. 


a  =  s-nu 

X.  =  <ja-r1-r2-...-r._1>, 
x  =  min(n,  t-r,-r7-. . ,-r._?). 


i-1. 


j  =  1,2,..., 
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2  n 

(ii)  g^t.u)  *  l  (c)Qt(t-s»s,u) 
s-s^ 

where  s1  and  s2  are  defined  in  (ii)  of  Theorem  2.2  by  taking 
y  =  nu. 

Theorem  2.3.  (i)  If,  for  given  P*(0  <  P*  <  £^y), 

g^t^.u^)  >_  P*(n^+n)  observing  that  t  =  ( t^  ,t2»  -  -  -  »tk) »  then, 

«3(v|t)  satisfies  the  P*-condition,  where  v^_^+1  =  u..,  i  =  l,2,...,k. 
(ii)  If,  for  some  given  P*(0  <  P*  <  £yf).  g^tt- ,u^)  >,  P*(  ^  )  for  all 

ti  =  0,1 ,2,. . . ,(i+l)n  and  i  =  1,2,. ..,k,  then  62(u* )  satisfies 

the  P*-condition  where  =  max  u. (t). 

k  1  0<t<(i+l)n  1 


Proof:  It  follows  from  (2.14),  (2.16),  (2.25)  and  the  definition  of 
63(y|t)  that  (i)  holds.  For  (ii),  we  note  that  P(A(i;  tj.t^)) 

(defined  by  (2.16))  is  increasing  in  and  (2.26)  holds  for  all  t1  *  0 
...,(i+l)n.  Finally,  it  follows  from  (2.14),  (2.18)  and  (2.2.5). 
Finally,  we  define 


S’.(t;ara2;B,Y)  =  {(x]  .Xg,. . .  ,Xf  )|a]  <  X}  <  c*2,  o^+B  <  x]+x2  <  o2+y. 

i-1  i 

. . .  ,o-|+(i-2)6  1  ^  xj  io2+(i-2)Y»  ^  Xj  *  t) 

c;(t;y,«)  *  {(x0,x1,...,x1)|x0-Y  <  x1  <  x0+6,2(x0-Y)  <  Vx2  < 

i  i 

2(x0+d),...,i(x0-r)  <  Jx.  <  i(xQ+6),  J  xf  -  t> 


<;  ><?  )...(!  ) 

Cj(t;a-|,a2;B.Y)  1  2  1 


uj(t;ai,o2;6,Y)  *  l  _  . x 


vi(t;Y.«)  -  i  (nx)(nx )...(; ). 

1  cj(tiY,6)  X0  X1  X1 
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By  the  analogous  arguments,  we  obtain  the  following 

[<*2]  s£ 

Corollary  2.4:  (i)  ujU'.a.  .cu'.e  ,y)  -  l  ("){£("  ){••• 

r|"<(X|^  1  r2~s2  2 


{  IJ  ("  )...{  I  (J  )L  r  "  _  )}. 

_  r .  _  r.  ,  't-r,-r,~. .  ,-r.  , 

rj“sj  J  ri-i'si-l  1-1  1  2  1_1 


where 


Sj  -  ^  ”  2,3,. • .  ,1-1. 

sj  =  [a2+(j-l)Y-r1-r2-...-rJ._1] 
x  =  min{n,  t-r,-r2-. . .-r^_2} 

Sn 


’1 


(ii)  v-(t;y,6)  =  l  (J)uJ(t-s;s-Y,s+6;s-Y,s+6) 
s=s0 

where 


sQ  *  maxtC^Jy*-],  t-in,0> 
s2  =  min{<|^->,  n,t}. 


Computations  of  ( t^ )  and  u! 

For  given  n,  i,  P*  and  t.  (0  <  t-  <_  n(i+l)),  we  may  start  by 
taking  u  =  0,  then  apply  Corollary  2.3  to  compute  g^(t^,0).  If 

g^t^O)  >_  P*(n^+n),  then  stop  and  take  u^(t^)  =  0.  Otherwise,  increase 

u  by  an  amount  of  y^-,  l.e.  u  =  jpj-.  Again,  using  recurrence  relation, 
compute  g^t^.u).  If  g^(t.,u)  >  P*(^+n),  then  we  take  u^t^)  =  u. 

Otherwise,  we  continue  this  process  until  for  the  first  time 
gi(t1,u)  >  P*(1J^n).  Then,  we  take  Vi+l^k-i+l*  "  “i^k-l+l^ 


L. 


1 


'  rlttifflMMyBliailliiii  -  .  I  MMMjjWnafri 
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For  the  procedure  6,,  we  take  ul  *  max  u^t*).  Some  special 
1  1  0<t.<(i+l)n  1  1 

values  of  u^(t^)  associated  with  n(=  5(1)10),  P*(=0.90,0.95)  and 
i(=  1(1)5)  are  tabulated  in  Table  2. 


3.  Some  Comments  on  the  Conditional  Isotonic  Rule  6^ 

As  it  can  be  seen,  the  unconditional  selection  procedure 
defined  in  section  2.2  always  satisfies  the  P*-condition.  As  a  matter 
of  fact,  it  follows  from  (2.14)  that  the  infimum  of  the  probability 
of  a  correct  selection  attains  1  if  u.  =  1.  On  the  other  hand,  when  the 
condition  of  total  sum  of  observations  is  imposed  on  the  event  B(i;  u.) 
defined  by  (2.17),  it  becomes  the  event  A(i;  t^,u^)  defined  by  (2.16), 
and  its  probability  of  A(i;  t^,u.)  decreases.  It  follows  from  (2.18) 
and  (2.25)  that 


(3.1) 


n(i+l)  g^.u)  i 

=  L 

j*0 


t.=0 


TnT+rT 
'  j 


vi  ■ 


where  Yq,Y^,...,Y^  are  iid  such  that  Yq  is  b(n;p)  for  some  unknown 
0  <  p  <  1.  Let 


fm+n> 
j 


(3.2)  hi(j;  u)  =  g^J.uJ/r1*"),  j  «  0,1 ,2,. . .  ,(i+l)n 


where 


0  <  u  <  1. 
Then,  (3.1)  becomes 


(i+l)n 

(3.3)  P(B(i ;  u))  *  l  h.(j;  u)«b(j;  (i+l)n,p) 

j«0 

i 

where  b(j;  n;  p)  is  the  probability  of  the  event  Yr  *  j).  If 


i .  .  i  . 


h.j(j;  u)  >  P*  for  all  j  =  0,1 ,2,. ..  ,(1+l)n  and  for  some  0  <  u  <  1,  then 
P(B(i;  u))  >  P*.  However,  it  is  not  true  that  h^(j;  u)  >_  P*  for  all  j, 
when  the  right  hand  side  of  (3.3)  is  not  less  than  P*.  As  some  computa¬ 
tions  show,  for  some  j  (when  j  Is  large),  h..(j;  u)  never  reaches  P* 

(e.g.  P*  a  0.95)  no  matter  how  large  u  is  .  This  undesirable  situation 
fortunately,  never  occurs  in  the  conditional  selection  procedure 
proposed  in  [3]. 


Table  1  devalues  (i  *  1(1)5) 

In  this  table,  the  smallest  d.-values  (i  -  1(1)5)  satisfying 

ai (**1 )  -  P*  defined  *>*  Lemma  2,1  are  tabu1ated  for  Po  *  O.KO.no.5. 

For  values  of  pQ  *  (0.6(0. 1)0. 9) ,  the  problem  can  be  treated  by 
considering  failures  instead  of  successes.  In  the  table,  the  upper 
entry  is  associated  with  P*  *  0.90  and  the  lower  associated  with 
P*  =  0.95.  The  entries  under  the  column  of  (O.a  -  O.b)  mean  that 
the  devalues  are  the  same  for  all  these  pQ-values  in  the  range 
(O.a  -O.b).  A  "  -  "  means  the  same  value  as  the  preceding  value 
in  the  same  column.  A  "  *  "  means  the  same  value  as  the  preceding 


value  in  the  same  row. 


Table  1 

d, -values 

devalues 

devalues  ** 

d.- values 

0.050  0.200 

0.045  0.227 

0.042  0.167 

0.250 


•A  missing  entry  In  this  table 
means  that  the  value  of  dj  Is  tl 

same  as  the  associated  pQ.  And 
also  the  values  of  dj  In  the  column 
of  pQ  «  0.1  areall  0.1. 


0. 1-0.4 

0.5 

0.100 

0.300 

0.083 

0.417 

0.071 

0.214 

0.357 

0.063 

0.250 

0.056 

0.278 

0.050 

0.200 

0.300 

0.045 

0.227 

0.042 

0.250 

0.038 

0.192 

0.269 

0.036 

0.214 

0.033 

0.167 

0.233 

0.031 

0.188 

0.219 

0.029 

0.176 

0.206 

0.028 

0.167 

0.222 

0.026 

0.184 

0.025 

0.150 

0.200 

0. 1-0.4  0.4  0.5 

0.067 

0.067 

0.133 

0.333 

0.056 

0.111 

0.389 

0.048 

0.333 

0.381 

0.042 

0.250 

0.333 

0.037 

0.296 

0.033 

0.200 

0.300 

0.030 

0.242 

0.028 

0.222 

0.250 

0.026 

0.205 

0.282 

0.024 

0.214 

0.022 

0.200 

0.244 

0.021 

0.188 

0.229 

0.020 

0.176 

0.216 

0.0)9 

0.185 

0.222 

0.018 

0.193 

0.017 

0.150 

0.200 

**A  missing  entry  In 
this  table  In  column 
of  p0  •  0.4  means  the 

value  of  dj  Is  the 

same  as  the  associated 
value  In  column 
0. 1-0.4. 


DEB 

0.4 

0.5 

0.050 

0.100 

0.150 

0.300 

0.400 

0.042 

0.125 

0.375 

0.036 

0.107 

0.321 

0.393 

0.031 

0.094 

0.281 

0.313 

0.028 

0.083 

0.278 

0.333 

0.025 

0.075 

0.300 

0.325 

0.023 

0.682 

0.250 

0.318 

0.021 

0.833 

0.250 

0.292 

0.019 

0.577 

0.250 

0.269 

0.018 

0.536 

0.267 

0.286 

0.017 

0.050 

0.233 

0.250 

0.016 

0.047 

0.234 

0.250 

0.015 

0.059 

0.235 

0.250 

0.014 

0.042 

0.222 

0.250 

0.013 

0.053 

0.237 

0.250 

0.0)3 

0.038 

0.237 

0.250 

Table  2  (t^-values  (i  =  1(1)5) 

Associated  with  n  -  (5(1)10),  i'"*=  (1(1)5)  and  P*  *  (0.90,0.95) 

the  entries  are  u^(t^)-values  satisfying  g^(t^,u^)  >  P*(^^n),  defined 

i  1 

by  Corollary  2.3,  with  t^  =  ^  x j ,  the  total  sum  of  (i+1)  populations. 

The  upper  entry  is  for  P*  =  0.90  and  the  lower  for  P*  =  0.95.  A  "  *  " 
means  the  same  as  the  preceding  value  in  the  same  row  and  a  "  -  " 
means  the  same  value  as  the  preceding  value  in  the  same  column.  The 
constant  v.  =  v.(t.)  needed  for  the  rule  6,  is  given  by  v(t)  =  u.(t.)- 
It  should  be  pointed  out  that  missing  entries  in  the  table  were  either 
not  computed  or  have  been  dropped  for  sake  of  convenience. 
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Table  2 
N  =  5 


i/t 

l 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 

0.400 

0.600 

0.400 

0.600 

0.800 

o 

0.300 

* 

0.500 

★ 

* 

★ 

* 

C 

- 

0.500 

* 

0.600 

0.500 

★ 

0.700 

0 

0.267 

★ 

* 

★ 

* 

★ 

0.400 

0.467 

0.667 

J 

- 

0.467 

- 

- 

- 

- 

- 

- 

- 

4 

0.250 

★ 

0.150 

★ 

* 

0.250 

* 

* 

0.450 

0.850 

5 

0.240 

★ 

* 

★ 

* 

* 

* 

* 

* 

* 

N  =  6 

i/t 

1 

2 

3 

4 

5 

6 

7 

8 

9 

_ 10 _ 

1 

0.333 

0.500 

0.333 

0.500 

0.667 

0.500 

- 

- 

0.667 

- 

- 

- 

9 

0.250 

* 

0.417 

★ 

★ 

* 

0.500 

0.583 

* 

C 

- 

0.417 

- 

0.500 

0.583 

* 

0.583 

★ 

0.833 

0.222 

* 

0.333 

0.278 

0.222 

★ 

0.333 

0.278 

0.389 

* 

O 

- 

0.389 

* 

- 

- 

- 

- 

0.389 

- 

- 

yl 

0.208 

★ 

0.125 

* 

* 

★ 

0.208 

* 

* 

* 

- 

* 

- 

- 

- 

- 

- 

- 

0.375 

c 

0.200 

* 

* 

* 

* 

★ 

* 

* 

* 

* 

D 

- 

- 

0.367 

- 

- 

- 

- 

- 

- 

- 

Is*  o 
in  o 

o  co 
cr>  co 

cn 

rs. 

co  in 

do 

—  co 

d  d 

♦  r— 

o 

fs^ 

o> 

In 

w  i 

d 

co 

CM  % 

d 

o 

r~  1 

d 

00 

QO 

—  »o 

N  CM  O 

m  i  m 

8 
*  CM 

f". 

.  2 

M 

Z  P> 

o  o  o 

d 

d 

0.429 

* 

* 

«  1 

8s: 

o  o 

dd 

CO 

0.571 

* 

0.500 

061‘0 

CV 

o. 

o 

in 

0.429 

* 

0.429 

0.238 

* 

*  » 

0.286 

0.571 

0.357 

0.286 

* 

0.107 

co 

0.429 

* 

0.357 

* 

0.333 

«  • 

CM 

1.286 

1.214 

1.190 

1.179 

« 

oo 

dd 

*  CM 

o 

CM 

Tab 
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